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Abstract 

We show an explicit attempt to interpret the multi-muon anomaly recently claimed by 
the CDF collaboration in terms of a light scalar singlet <p which communicates with the 
standard quarks either through a heavy scalar or a heavy fermion exchange. Building on 
arXiv:0810.5730, that suggested a singlet 4> with a chain decay into a final state made of four 
rf pairs, we can simulate most of the muon properties of the selected sample of events. Some 
of these properties adhere rather well to the already published data; others should allow a 
decisive test of the proposed interpretation. Assuming that the test is positively passed, 
we show how the PAMELA excess can be fitted by the annihilation of a TeV Dark Matter 
particle that communicates with the Standard Model via the new light singlet(s). 



1 Introduction and statement of the problem 



The CDF Collaboration has recently published a study of a significant sample of multi-muon 
events with unexpected properties pQ. The DO Collaboration performed a similar search without 
finding similar events [2]. While it is still possible that these events will be in the end accounted 
for in terms of standard physics and detector effects, yet a small but significant fraction of them 
has characteristics which are peculiar enough to deserve attention. This has in fact prompted 
a phenomenological conjecture [3] that tries to explain them in terms of new physics: the pair 
production of a relatively light Standard Model singlet, 0, each with a cascade decay into a final 
state made of 8 tau leptons. An intriguing feature of this explanation is that it describes well the 
sign-coded multiplicity distribution of the additional muons, from two on, found in a cone^] around 
the direction of a primary muon. In this interpretation, such a distribution is simply related, given 
the detector efficiencies in muon identification, to the r — > // decay probability. On the other hand, 
the apparent difficulty in exhibiting a plausible production mechanism of the 0-pair casts doubts 
on the significance of the proposal and, even more importantly, prevents a decisive comparison of 
the model with the data. In spite of its uncertainties and ambiguities, this situation motivates us 
to try to go further. 

A main problem one faces right at the beginning is how to select, from the overwhelming 
number of background events, those ones that may represent the signal. Secondly, a proper 
understanding of the event properties that involve hadronic tracks requires a detection simulation 
that we cannot do. To get around these problems we adopt the following strategy. We first 
concentrate our attention on a sub-class of about 4000 events that contain at least two muons 
(i.e. one primary and at least one additional) in each of the two cones. These events have also 
been highlighted in [Tj [3] and found to have several properties especially difficult to understand 
in terms of known physics. We assume that these events constitute the signal and we seek for a 
production mechanism of the 0-pair that can account for the measured properties of the muons 
contained in them. Out of all possibilities, discussed in Section 2, somewhat surprisingly one single 
effective operator of dimension 5 fits the invariant mass distribution of all muons in this sample 
for a definite value of its scale. The operator is of the form qq<f) 2 , where q is a first generation 
quark. 

Encouraged by this result, which fixes the production mechanism, we analyze in Section 3 
a few other measured features of the multi-muon events: the ^ multiplicity, impact parameter 
and invariant mass of muons in one cone. Given our definition of the signal, the background can 
be clearly identified. We show that it should be possible, with the data at hand, to fit the tail 
of the muon impact parameter distribution by adopting the chain decay model proposed in [3], 
<p — > 20i — > 402, where 0i,2 are two other SM singlets and 02 decays into a tau pair with a lifetime 
around 30 picoseconds. Our main purpose is to make possible a close comparison of the model 
with the data so as to allow a decisive test. To this end we discuss in Section 4 a few quantitative 
predictions of the model that should be compared with the data. 

1 The events in question are triggered by a pair of muons with px > 3 GeV and \rf\ < 0.7. Additional muons 
with pt > 2 GeV and |?7| < 1.1, which are often present in these events, are grouped into two cones of opening 
angle 36.8° around the directions of the trigger (also called primary) muons. Muons which happen to lie outside 
of these cones are ignored. 
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The 0-pair production cross section needed to explain the signal is about 200 pb, in turn 
requiring a scale of the effective operator that describes it of about 100 GeV. To make it acceptable, 
we find it necessary to deconstruct the effective operator in terms of renormalizable interactions, 
with new particles mediating the communication between the quarks in the proton and the hidden 
sector that contains 0. These new particles and interactions must not lead to any unseen structure 
in the multi-muon data and must be consistent with known experimental constraints. Given the 
size of the cross section and the low effective scale involved, this proves nontrivial to achieve. 
Nevertheless we illustrate in Section 5 two minimal models that can pass the test, to the best that 
we can tell. One involves a scalar exchange in the s-channel and another a fermion exchanged in 
the t-channel. In both cases the couplings to the singlet saturate perturbation theory. We also 
briefly comment on the problems that one may face in turning them into complete satisfactory 
extensions of the SM. 

If this interpretation of the CDF events will resist a further scrutiny along the lines we are 
proposing, the existence of the light scalar sector suggests a connection with recent astro-particle 
data. In Section 6 we consider a hidden-sector Dark Matter model where DM annihilations into 
0's that decay into r's can provide the positron excess claimed in e + cosmic rays by PAMELA [4], 
while giving in the e + + e~ spectrum a feature somewhat smoother than the peak claimed by 
ATIC [5J. Conclusions are summarized in Section 7. 

2 The production of the signal events 

As mentioned, we first concentrate our attention on the events with two cones containing at least 
two muons each, which we call signal events as their features are peculiar enough that they should 
be minimally polluted by known backgrounds and could therefore be a quasi-pure sample of a 
beyond the SM signal. For an integrated luminosity of 2.1 fb _1 there are about four thousand 
such events pQ. If they have to arise from — > 4tt, this requires a significant pp — > 00 cross 
section, above 100 pb. Furthermore the invariant mass distribution of all muons contained in these 
events, Fig. 35a of pQ, which we aim to explain, does not show any special feature other than a 
threshold rise and an extended tail. In view of this we are led to consider an effective operator 
bilinear in and with two gluons or a quark-antiquark pair. There are three such operators of 
dimension less or equal to (rj 



where q is a quark field, either u or d, and G the gluon field. We will normalize our plots below 
assuming that O5 couples only to u, while Oqf couples to both u and d with equal strength. 

With the amplitude corresponding to each of these operators, we have simulated the 0-pair 
production followed by the decay chain 




(2.1) 



0^2 0! 



4 2 -> 8 r 



(2.2) 



2 For definiteness, we assume that <p is a complex field. 
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with = 15 GeV; see [3] and Section[3]for a discussion of this choice. Imposing the experimental 
cuts and detection efficiencies^] we get for each of the above operators the signal cross section^ 



n ,200 GeVV , / 200 GeV \ 4 , / 200 GeV Y 

a 5 = 0.38pb - , a 6F = 1.3pb , a 6G = 0.52pb' ^ 



— J , (7 6F = 1.3pbl J , cr 6G = U.52pb ^ J 

For reference the total pp —>■ 00* cross sections (without imposing cuts and efficiencies) are 
<x 5jtot = 220pb, cr 6Fitot = 105 pb, cr 6G , tot = 370 pb for A = 200 GeV. 



Fig. [T] shows how well each of these three production models can account for the experimental 
distribution of the invariant mass of all muons contained in the two signal cones, shown in Fig. 35a 
of [1] . The hardness of the observed distribution motivated us to consider production via effective 
non-renormalizable operators. The single relevant parameter here is the effective scale attached 
to each individual operator. We find it remarkable that a single operator, O5, can fit the full 
distribution. For definiteness we take q = u only, which requires A = 85 GeV (fixed in this case to 
reproduce the total number of events). On the other hand Oqf and Oqg can separately reproduce 
only the tail and the low-mass region close to threshold, respectively. 

The dashed red line in Fig. [IJi shows the spectrum shape which results if 00* are decay products 
of a narrow 300 GeV scalar resonance: qq — > S — > 00* (in arbitrary normalization). The purpose 
of including this spectrum, clearly unable to reproduce the experimental data, is twofold. First, 
this case was considered in and our simulations agree. Second, since the 300 GeV ^-exchange 
gives a spectrum peaked around 30 GeV, we can deduce that muons carry on average around 1/10 
of the total 00* invariant mass. In particular, events at the tail of the experimental spectrum, 
M muons ~ 100 GeV, should correspond to V§ ~ 1 TeV. 



3 Some measured properties of the multi-muons 

Concentrating on the 0-production described by O5 with A fixed at 85 GeV, we here try to 
reproduce a few other significant properties of the muons in the CDF events of [I]: 1) the already 
mentioned sign-coded multiplicity distribution; 2) the invariant mass of the muons in a single 
cone; 3) the muon impact parameter distribution. Notice that data of [lj about 1) and 3) also 
contain events outside the sub-sample discussed in Section [2] that we view as signal events; thereby 
we allow for the presence of a background. 

3 Trigger muons must have transverse momentum pr > 3 GeV, rapidity \r)\ < 0.7, and are detected with efficiency 
PcS = 0.44; additional muons must have px > 2GeV, \t)\ < 1.1 and p c g — 0.838. Furthermore, the invariant mass 
of the trigger muon pair must be between 5 and 80 GeV; the relative azimuthal angle for the opposite sign trigger 
muon pair must be below 3.135 radians. 

4 Our simulations were performed by adding the necessary new particles and processes in Pythia 8.108 [B] and 
were cross-checked with the help of home-grown Monte-Carlo programs written in Mathematica (both stand-alone 
and passing Les Houches events to Pythia for parton shower). 
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Figure 1: The invariant mass of all muons in events with at least two muons in both 
cones. The experimental data are from Fig. 35a of [1], L int = 2.1 fb^ 1 . 

3.1 Sign-coded muon multiplicity 

For reasons of graphical illustration, the sign-coded muon multiplicity M. is defined for each of 
the two 36.8° cones around the trigger muons by the formula 

M = N os + 10N SS 

where Nqs {Nss) is the number of additional muons in the cone having opposite (same) sign as 
the trigger muon. Thus e.g. M. = corresponds to cones without additional muons. In Table 1 
and Fig. [2] we compare the data with our simulation. 

The shape of this distribution is largely independent of the production mechanism, and indeed 
our Fig. [2] is quite similar to Fig. 1 of [3] . Since we have fixed the production mechanism (the O5 
operator), we can simulate not only the shape but also the absolute normalization (A = 85 GeV), 
while in Fig. 1 of [3] normalization has been fixed arbitrarily. We see in particular that we can 
reproduce essentially all events (within errors) with 3 or more muons in a cone: only a relatively 
small systematic excess of not understood backgrounds is here possibly needed to fully account 
for the event rates. However, only about a third of events with 2 muons in a cone is reproduced. 
Since we fixed the normalization by using the events with two or more muons in both cones, this 
shows that events with only 2 muons in one cone and 1 muon in the other are still significantly 
contaminated by backgrounds. 

The fraction of events without additional muons that we are able to reproduce is negligible. 
Almost all of these events must be background. Indeed according to pQ at least 50% of events 
without additional muons can be due to in-flight decays of kaons and hyperons producing real 
muons or pion punchthroughs. 
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Figure 2: Comparison of our production model with the experimental data (Fig. 22b of 
significant part only). Every cone enters separately into this histogram (thus every 
event enters twice). 
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Table 1: Numerical content of Fig. Experimental data taken from Table X of ^ 
(significant part only). The data correspond to L int = 1426 pb^ 1 . The statistical errors 
of our simulation f~ 0.2y/Nbi n ) are negligible compared to the experimental errors. 
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Nt,IN>2>2 as a function of p 



Na/N>2>2 as a function of p 
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Figure 3: The numbers of events with 3 muons in a cone N 3 (the sum of bins 02,11,20 
in Table 1) and with 4 muons in a cone iV 4 (the sum of bins 03,12,21,30) relative to the 
number of events with at least two muons in each cone N>2,>2, simulated as a function 
of the varying r — > \i branching ratio p. The bands show experimental values. The 
physical value p = 0.17 gives reasonable agreement in both plots. One can see that the 
ratio N 3 /N> 2 ,>2 is rather insensitive to the variation of p, while the ratio A^ 4 /A^> 2i >2 
prefers the values of ' p ~ 0.2, close to the physical value. 




Notice that the above agreement of simulation with experiment, in bins with > 3 muons, relies 
on the precise value of p = BR(r — > fx) m 0.17, thereby indirectly supporting the r interpretation. 
This can be seen by generating toy Monte-Carlos in which p varies in the 0.1 -j- 0.3 range (while 
keeping fixed the muon detection efficiencies). The number of events with 3 muons in a cone 
(including sign-coding), relative to the events with at least two muons in both cones, is largely 
independent of p. This is not unexpected, since the probability to put 2 additional muons in the 
same cone is roughly the same as to put them in different cones. On the contrary, the relative 
number of events with 4 or more muons in a cone was quite sensitive to p, to the extent that 
p = 0.1(0.3) led to clear shortage (excess) of events in these bins, compared to experiment, see 
Fig. |} At the same time the relative distribution of these events into the various sign-coded bins 
is determined by simple combinatorics and, as a result, proves to be insensitive to the variations 
of p. 



3.2 Single cone muon invariant mass 

The multi-muon events are characterized by low invariant mass of muons in a single cone, with 
a spectrum which sharply cuts off at 2 3 GeV. It was shown in j3] that this cutoff can be 
reproduced within the hypothesis of the decay chain (2.2) if the mass of (f) is not much above the 
8r threshold (14.2 GeV). This conclusion is largely independent of the production mechanism, 
and in particular holds within our model. In Fig. [4] we compare the data with our simulation for 
= 15 and 17 GeV (for events where both cones contain at least 2 muons, i.e. the signal events 
of Section [2]). 

The overall normalization is fixed in both cases to reproduce the total number of events. 
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Figure 4: Invariant mass , M cone , of all muons in a cone when both cones contain at 
least two muons. Data reproduced from Fig. J^la of 

While the data are reproduced in an acceptable way for these two masses, higher values (see [3] 
for rrifp = 20 GeV) start showing a clear deviation. 

We however see that a discrepancy is present at small M cone : while the data remain almost 
constant down to the lowest bin, the simulation drops there. The physical reason for this drop 
is intuitively clear: it corresponds to a small probability to have an almost collinear muon pair. 
Without a precise understanding of the detector (e.g. possible punch-through of hadrons mimicking 
muons) and of the statistical correlations between the various bins we cannot say more on this 
issue nor perform a precise y 2 analysis. 

It is legitimate to ask if the experimental distribution in Fig. [4] also contains a peak corre- 
sponding to the possible subdominant 2 decay channel, directly into muons, 2 — ¥ Such 
a peak would have to be located at m^, i.e. slightly above 2m T ~ 3.55 GeVj^J Indeed we see some 
excess of events around this value of M cone (which could become statistically significant with more 
statistics), but, at the same time, we do not forget an important source of dimuon pairs — cc 
mesons — present in the 3 -r 4 GeV mass region. From the model-building point of view, the 
branching ratio 4> 2 — * M + At~ ma y well be suppressed to an unobservable level (see Section [5]). 

3.3 The muon impact parameter 

For a given charged particle track, the impact parameter d is defined as the distance between the 
track and the primary interaction vertex in the transverse plane, see Fig. [5} 

One of the most striking properties of multi-muon events is that the impact parameter distri- 
bution of muon tracks has an exponential tail. Such a tail can be explained if the tracks originate 

5 And not at 7.2 GeV as in Fig. 1 of [7]. 
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Figure 5: The impact parameter d is the distance between the track projection on the 
transverse plane track pro j and the primary interaction vertex P, which is experimentally 
known with high precision (in particular due to many soft hadronic tracks originating 
from it), while the secondary vertex S is not measured unless more than one track 
originates from it. 

from decays of a long-lived particle. Notice that a boost 7 of the decaying particle does not affect 
the typical value of the impact parameter, since it increases the typical decay distance I ~ jctq 
but also decreases the typical collimation angle of decay products by the same factor 7. As a result 
the decay scale of the impact parameter distribution is of the order of the decaying particle's cr . 
The precise proportionality coefficient in general depends on the kinematic cuts. 

In Fig. [6] we compare one of the several published impact parameter distributions (impact 
parameters of primary muons of all cones containing exactly 2 muons) with our simulation. We 
assume that <fi — > 2(f) 1 — > 4 2 decays are prompt, while 2 — > rf has a long lifetime. As discussed 
in this assumption is forced on us by another piece of experimental data which we do not 
discuss here: no significant correlation between the impact parameters of different muonic tracks 
in the same cone. 

We see that both the decay rate and the normalization of the exponential tail can be reproduced 
for a 4> 2 lifetime around 30 psj^] Note that 'naive' decay exponent would give a 20 ps lifetime; 
the fact that the naive ultra-relativistic estimate underestimates the true lifetime by 30% in this 
particular case has already been noticed in [3] on the basis of a toy simulations. We also see that 
the large number of events at small d are not reproduced by our signal model. Most of these events 
must be due to QCD contributions, e.g. sequential b quark decays which had been subtracted in 
the distributions considered previously but not in Fig. [6j All data are extracted from PQ, so that 
we follow their definitions. 

4 Possibilities for further experimental tests 

Now that a concrete proposal for the 0-production via the O5 operator has emerged, what other 
possible tests could provide further checks of the model? Many such tests can be done using 
information about hadrons present in the multi-muon events. However, any study with hadrons 
would require detector simulation, especially because the relevant experimental information pQ 
involves hadronic tracks rather than, say, isolated jets. 

This value has to be taken with a grain of salt since the tail of the experimental distribution may be significantly 
contaminated by K s in-flight decays (lifetime 90 ps). 
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Figure 6: Primary muon impact pa- 
rameter distribution for cones con- 
taining exactly two muons. Data 
taken from Fig. 25a of /!]/. 



Figure 7: The simulated distribution 
for L xy for r = 30 ps and for dimuon 
secondary vertices (in arbitrary nor- 
malization). 



In this Section we discuss 3 further experimental tests which involve only muons and which 
might be relatively easier to implement. The corresponding experimental distributions are not yet 
available. These tests are: 1) dimuon displaced vertices; 2) muon track alignment; 3) deviation 
from back-to-backness in the CM frame of the hard process. For illustrative purposes only we will 
show in all cases the expected distributions from our simulation. 

All our predictions below are computed using our simplified simulation of the triggering and 
reconstruction process, see footnote 1. In order to perform a precise comparison with the data, 
the predictions should be recomputed using the full simulation of the detector, available only to 
the experimental collaborations. 



4.1 Dimuon displaced vertices 

If the model is right, a fraction of events must contain a displaced secondary vertex with two muon 
tracks originating from it. Such a vertex occurs when both taus in <p2 —> tt decay into a muon. 
According to our simulation, for Lj nt = 2.1 fb^ 1 there will be about 2000 such secondary vertices, 
for <p2 lifetime of 30 ps almost all of them within the CDF inner silicon vertex detector (10.6 
cm radius). Taking into account excellent tracking capabilities of the CDF detector, a significant 
fraction of these vertices should be identifiable. 

A useful experimental quantity to describe the secondary vertex distribution is called L xy , 
the distance between the secondary and primary event vertices projected onto the transverse 
momentum of the two-track system. In Fig. [7] we give our simulated distribution for L xy when 
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the secondary vertex is determined by two muon tracks. 



4.2 Muon track alignment 



As discussed in Section |3j cannot be much above 15 -f- 17 GeV. On the other hand we must 
have m^ 2 > 2m T otherwise the decay of 02 will proceed via off-shell taus with a typical lifetime 
much bigger than the needed 30 ps. One consequence of such restricted kinematics is that the 
four 2 are non-relativistic in the rest frame of the parent 0, with typical kinetic energy 

£kin = ^ - ~ (0.2 4- 0.7) GeV. 

As noticed in [7j , this leads to 02 decay vertices collinearly aligned along the momentum. Indeed, 
their typical separation S in the orthogonal direction is 

5~f3 -ctq, /?~(2E kin /m ( ^ 2 ) 1 / 2 ~0.3-0.6. (4.1) 

Since we are here concerned only with muon tracks, the 02 decay vertices cannot be recon- 
structed in most cases. However, a nontrivial experimental measure of their collinear alignment 
can be defined purely in terms of muonic tracks, for events with at least 3 muons in a cone. We 
will use the following quantity: 

. 1/2 

A = min ( Y^i d i( 

where di is the distance between the track of the z-th muon and a straight line passing through 
the primary vertex in the direction p. By tracks here we mean idealized reconstructed tracks, 
where the effects of bending in the magnetic field of the detector have been deconvoluted. Thus 
A = exactly when all tracks intersect a single line pointing to the primary vertex. Generically 



we expect A positive and of the order of Eq. (4.1), although this is likely to be an overestimate 



since the definition of A makes it automatically vanish for cones with only two muons. 

In Fig. [8] we plot the simulated distribution of A of events with at least 3 muons in a cone. 
A less steep distribution, similar to the one of the impact parameter d in Fig. |6j would instead 
arise from non-aligned r decays. In Fig. [8] we assumed two values of the mass: 15 and 17 GeV, 
777.02 = 3.6 GeV and = 7.3 GeV. 



4.3 Deviation from back-to-backness 

While the total muon momenta in each cone pip is only a fraction of their parent momentum, 
so that the muon invariant mass gives a poor information on the mass, pi is very strongly 
correlated to the direction of 0, allowing to test details of the production mechanism. 

The distribution of 0's in rapidity does not contain much information, being mostly shaped 
by the experimental cuts and by the parton distributions. More interesting is the distribution of 
the relative azimuthal angle Atp (i.e. the relative angle between p\ and P2 in the plane transverse 
to the beam), since the deviation from the back-to-back configuration (Aip = tt) can be due only 
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Figure 8: Simulated distributions of 
the muon track alignment parameter 
A, for events with exactly 3 muons in 
a cone. 



Figure 9: The expected 'deviation 
from back-to-backness' angle n — A(p 
distribution for events with at least 2 
muons in each cone. 



to QCD radiation in our model. In our model is a color singlet, so that only Initial State 
Radiation is present, leading to small deviations from back-to-backness. Fig. [9] shows our result 
for the simulated distribution of A<p for signal events containing at least two muons in each cone: 
the average angle is 

(180° - Ay?) » 7°. 

We relied on Pythia 8.108 with its standard settings of ISR parameters and no detector sim- 
ulation, and we have not attempted to quantify QCD uncertainties on this result. Once these 
issues are settled, such distribution can discriminate from other production mechanisms or pos- 
sible backgrounds. For instance, <ft production via the alternative gluon operator Oqg gives ~ 2 
times wider distribution since Initial State Radiation off gluons is stronger. A wider distribu- 
tion would also be obtained in presence of Final State Radiation, which could be the case if the 
multi-muon events are background associated with dijet events. Indeed, analogous experimentally 
measured distributions for bb production [1U] are significantly wider. 



5 Explicit realizations for O5 

We are not concerned here by the nature of the hidden world, and we do not pretend to explain 
the interaction that leads to the last step of the chain decay: 02 ~^ tt. Given its small width, we 
can attribute it to another effective operator 

^decay = H°(j) 2 (5.1) 
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where H s is the standard Higgs field, l 3 is the third-generation lepton doublet, r the singlet. The 
scale A' associated with it is much above the Fermi scale, unlike the A associated to 0-production. 



On the contrary, to be credible, the effective non-renormalizable operators in eq. (2.1) sup- 
pressed by a scale A as low as one hundred GeV must be "deconstructed" in terms of some explicit 
particle exchanges. What can mediate the pretty strong communication between the SM particles 
and the hidden world represented by <fi or the other light particles to which </> decays, without 
running in conflict with known experimental constraints? 

We concentrate our attention to O5, leaving the two dimension-6 operators to Appendix A. 

There are only two ways in which the operator O5 can be generated by the mediation of a 
single particle: 1) a scalar exchanged in the s-channel, with the quantum numbers of the standard 
Higgs doublet, but different from it and with a negligibly small vacuum expectation value; 2) a 
colored fermion exchanged in the t-channel. We analyze both cases in turn. 

5.1 Scalar exchange 

The renormalizable Lagrangian that gives rise to O5 after integrating out the heavy scalar doublet 
H, with hypercharge —1/2, is 

L H = \ q quH + \H ] H s \(f)\ 2 - M 2 H H ] H , (5.2) 

where q is the first generation left-handed quark doublet and u is the right-handed up quark. 
(Similar considerations would hold for the down quark.) In this way, at energies below Mh, one 
generates an effective Lagrangian with O5 for the up quark and a scale 

which, to match with Fig. [TJ should be about 100 GeV. 

Requiring as the perturbativity limit that the partial width H — > cfxfi* be less than the H 
mass^j we get 

T{H - #*) ~ < M H Xv<Vl6^M H . (5.4) 

From the last two equations we conclude that 

M H /X q < 700 GeV. (5.5) 



This bound applied to the Lagrangian (5.2) raises issues of potential conflicts with experiments, 
which we briefly address in Section 5.3. Here we focus on its possible effects at the Tevatron itself. 
The first of these effects concerns the same pp — ► <fxp* production, since the effective-operator 



approximation in describing the if -exchange may be invalid. This is illustrated in Fig. 10 where 
the effective-operator result of Fig. [IJi is compared with the exchange of the if-scalar with a mass 
of 1 TeV and a width of 600 GeV. For these values of the parameters we see no conflict with the 
data, which may even be better described in the second case. 



7 For Mh < 3 TeV this turns out to be stronger than the NDA bound A < 167r 2 
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Figure 10: The analogue of Fig. where the contact term interaction generated by the 
local operator O5 is compared with a scalar s-channel effect and a fermion t-channel 
exchange generating the same operator at low energies. The couplings are slightly ad- 
justed so that the total number of events is the same in all cases. 



Another effect of the if -exchange might have shown up in dijet properties. CDF has published 
studies both of the angular distribution of the dijets [8] and of a possible resonant structure in 
the dijet mass spectrum [9]. In both cases it is nontrivial for us to compare these results with the 
signal we might expect. For the case of the angular distribution, the if -exchange generates the 
operator 

whereas the analysis in [8] concerns the operator (including the sign) 

1 ^2 



2A 2 



{QL'fQLf, Q = (u, d) 1 , A > 700 GeV at 95% C.L. . (5.7) 



The analysis leading to this bound is based on the fact that operator (5.7) interferes constructively 



with the QCD qq — > qq scattering amplitude, leading to an increase of events at large scattering 



angles. Our preliminary analysis shows that the interference term between the operator (5.6) and 



the QCD amplitude is smaller. Thus we expect that the bound in (5.5) should not pose a problem. 
All this is preliminary, however, also because the resonance may be within the accessible range. 
Indeed the same resonance search in [9] could be of importance to us. The problem in this case, 
however, is that the published limits concern only narrow resonances, when the width is within 
the resolution, while we are dealing with the opposite case. How to rescale from this, if possible 
at all, is to the least nontrivial. We believe however that the current results may only constrain 
Mh weakly enough to allow consistency with our effects in a significant range of the parameters. 
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Figure 11: Three ways to generate O5 via new heavy quark exchange in the t-channel. 



Only the last way, corresponding to the Lagrangian (5.9), is phenomenologically ac- 
ceptable. 



5.2 Fermion exchange 

The other way to generate the operator O5 is by a heavy quark exchange in the t-channel, see 



Fig. 11 In principle one can do it with one SU(2)x-singlet U or with one doublet Q, but this 



requires large couplings of the form 

X q qUH s , X q QuH s , (5.8) 

respectively for the U or the Q cases (here q is the standard SM left-handed doublet). 

However, after diagonalization of the full mass matrix in the 2/3-charge sector, this would 
lead to a large wrong component of the physical or w^ h , which is not acceptable. One needs 
therefore to introduce both a U and a Q and the Lagrangian 

L F = (X g qQ<P + X u Uu<f)* + XQ R U L H S + h.c.) + M Q QQ + M V UU, (5.9) 

which, in the local limit and the approximations Mjj m Mq >A«, gives rise to O5 with 

XqXuAV 

As in the case of a scalar exchange there are effects of non-locality, depending on the value of 
Mjj ~ Mq. The t term in the t-channel fermion propagator, (t — M^) _1 , while negligible for small 
t, suppresses the rate for events with the large invariant mass. The strength of this suppression 
can be estimated from t ~ — s/2 (relevant for a ~ 90° scattering angle). As discussed in Section 
|2j s ~ 1 TeV for events near the tail of the M muons distribution in Fig. [TJ Quantitatively, this 



effect is illustrated in Fig. 10 for M\j = 1 TeV, showing that a much lower value of the heavy 
fermion mass would lead to a too strong depletion of the tail in the invariant mass distribution 
of the muons. To comply with A w 100 GeV, Mjj > 1 TeV requires therefore X q X u X > 60, i.e. a 
collection of pretty strong couplings. 



The Lagrangian (5.9) generates also several 4-quark interactions at one loop, potentially rele- 
vant for the dijet physics. The typical scale of these operators will be ~ AttA w 1.3 TeV, which, 
as previously discussed, should be compatible with current limits. 
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5.3 Problems for a more complete theory 

The formulation of a full proper extension of the SM that incorporates in a consistent way both the 
hidden and the mediation sectors, in one of the two versions above, would obviously be premature 
at this stage and, as such, is outside the scope of this work. Here we briefly mention some of the 
problem one would have to face. They are of at least three kinds: 

1. Flavor problems; 

2. Naturalness problems associated with the lightness of (and the other hidden particles); 

3. In the scalar-exchange case, a naturalness problem associated with the need of a vanishingly 
small vacuum expectation value of H. 



The flavor problems are associated with the presence of the effective 4-quark interactions 
already discussed in connection with possible dijet signals. Even in the 2/3-charge sector the 
limits on the effective scale of these interactions, if not properly aligned in flavor space, are 
severe. One at least technically acceptable way out may consist in assuming that all the Yukawa 
interactions in Lh, Eq. (5.2), or Lp, Eq. (5.9), and in the SM be flavor diagonal except for the 
usual Yukawa coupling of the down quarks, which would therefore be the only source of flavor 
breaking. 

The issue of the naturally small masses of the hidden sector scalars has a lot to do with their 
nature and their possible internal structure, about which little is known. Nevertheless it may be 
in any case nontrivial to accommodate in a natural picture their large Yukawa couplings as in 
Q or in Q. 

Finally, from (5.2) a non vanishing vacuum expectation value of H gives rise to contribution 
to the up quark mass, 8m u = X q (H), which needs to be kept under control. In turn, since X q 
cannot be too small, this is a severe limit on the mixing squared mass term m 2 H if H s in the scalar 
potential. Seen another way, the operator O5 leads to a radiative contribution to the up quark 
mass from a loop. This radiative correction could be removed, without affecting the multi muon 
signal, by replacing |0| 2 — > <fi 2 in 5 and making proper adjustments in the deconstructed versions. 



6 Dark matter connection? 

Let us assume that the proposed interpretation of the CDF events is correct: a hidden sector 
with particles of mass of order 10 GeV, that decay into r's and coupled to the Standard Model 
trough new states Q, U or H with mass in the TeV region. Here we suppose that this TeV sector is 
augmented by a stable neutral state \ that is the dark matter. To connect to the multi-muon data, 
we also assume that \ is also strongly coupled to the hidden sector by the Yukawa interactions 
Xx(0, 0i, 02) for fermionic x, or the quartic couplings X*x(0*0, 0*0i, 0202) f° r scalar x, so that 
the dominant dark matter annihilation channels are 

XX - 0*0, 0i0i, 0^02 - (rr) 8 , (r r) 4 , (fr) 2 . (6.1) 
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Figure 12: The cosmic ray data about the positron fraction e + /(e + + e~) (left), e + + 
e~ (middle), anti-proton fraction p/p (right) compared to the expected astrophysical 
backgrounds (lower shaded curve) plus the Dark Matter excess (upper red curve). 

If this annihilation cross section is large enough, the e + from r decay are able to explain the 
excess positron cosmic ray signal reported by PAMELA [4], for m x > (2, 1, 0.5) TeV, for dominant 
annihilations to 0, 0i,02 respectively. The predicted spectrum for the positron excess is shown 



in Fig. 12, that applies equally well to annihilation via or 0i or 02 for m x = (4,2,1) TeV, 
respectively. Since the kinematics is such that the r are produced non-relativistically in the rest 
frame of the decaying hidden sector scalar, the positron spectrum is similar in the three cases. A 
smoother spectrum arises if instead DM has more than one annihilation channel with comparable 
branching ratios. (See [TT] and references therein for a discussion of the expected astrophysical 
backgrounds and uncertainties). 

The connection with the CDF multi-muon data is that the signals involve r decays. As a 
consequence there is no associated p excess, compatibly with PAMELA observations jl]. 

Future experiments measuring the positron excess at energies beyond 100 GeV should find 



that the excess persists smoothly to higher energies, with a turnover behavior, shown in Fig. 12 
that is not as sharp as when annihilations occur directly to electrons or muons. 

To obtain the PAMELA positron excess, the XX annihilation cross section must be larger than 
that required for a thermal freeze-out abundance of \ by a boost factor 

B « 10 3 ^, (6.2) 

independent of whether the annihilation yields 2,4,or 8 fr pairs. A contribution to this boost 
factor arises from a Sommerfeld enhancement factor, from a ladder of exchanges between the 
initial annihilating \ states, and is significant since -C m x . If x is a fermion the Yukawa 
coupling xx4> leads to this ladder, while for \ a scalar a new trilinear scalar coupling must be 
added, xx4>- F° r the case °f the Yukawa coupling, y%{S + il§A)xi the desired thermal freeze- 
out abundance results if y 2 ~ m x /1.5TeV. For m x in the range of 1 to 10 TeV the Sommerfeld 
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enhancement factor is approximately given by 



Tfl 

i?^ 200._2L (6 . 3) 



Comparing with (6.2), this is only a factor 5 below the required boost factor. Hence if local 
clump iness of the halo provide the extra ~ 5 enhancement, our theory is able to yield the PAMELA 
positron excess even with dark matter produced thermally. Such a local dumpiness is helpful in 
reducing tensions with limits on gamma [TT] and neutrino fluxes. 



In Fig. 12 we fixed the DM mass m x trying to reproduce the peak suggested by the ATIC [S] 
observations of the e + + e~ spectrum below 1 TeV. However, the spectral feature produced via 
r decays (as suggested by the CDF anomaly) is less pronounced than that arising from DM 
annihilations to electrons or muons, so that higher statistics data will test the hypothesis of 
annihilations via r pairs. ^ 



7 Conclusions 

We explored the possibilities for a new-physics interpretation of the CDF multi-muon events, 
adopting the phenomenological conjecture of |3] which links the muons to the pair production of 
two hidden sector scalars, each decaying into 8 r leptons. Within this conjecture we addressed the 
major puzzle left unexplained so far: a plausible production mechanism which could explain the 
total muon invariant mass spectrum. The hardness of this spectrum suggests to explore higher- 
dimension local operators connecting two quarks or gluons to a pair of 0's, and our systematic 
search revealed that one such operator, namely the dimension-5 5 = qq\4>\ 2 /A, reproduces the 
measured distribution for A ~ 100 GeV. 

Encouraged by this discovery, and by the ability to reproduce several other experimental 
distributions, we proceeded to discuss the next crucial question: is it at all conceivable that such a 
low scale can be generated in a more complete theory without producing any other unseen effect? 
We presented the two minimal ways of solving this problem, either by exchanging an s-channel 
new heavy scalar, or by exchanging a t-channel new heavy quark. In both cases 4-quark operators 
are unavoidably generated along with O5, with a scale of the order a TeV. While the comparison 
with the existing Tevatron limits from dijet studies is nontrivial, our preliminary analysis shows 
that in both cases no obvious contradiction exists. Furthermore we found that the effects due to 
momentum dependence of the exchanged particle propagators may not disturb the agreement of 
the effective-operator limit with the total invariant mass distribution. 

We think that all this offers a definite and consistent framework which can be subject to 
stringent tests by further comparison with the data. In fact, since the production mechanism is 
now largely fixed, it should be possible for the experimentalists to pronounce a final word on this 
model. We propose possible tests involving muons only; more tests involving hadronic tracks can 
be imagined. 



3 Note Added. The FERMI experiment [T5] did not confirm the ATIC peak and observed an e + + e spectrum 



compatible with the model prediction in figure 12 3 
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If we assume that the proposed interpretation of the CDF events is correct, further work in 
many different directions can and must be done. Here we have not resisted from drawing a possible 
connection between the CDF data and the putative Dark Matter signals seen in PAMELA and/or 
ATIC. Assuming that Dark Matter is a hidden sector particle that annihilates into the 0-scalars 
decaying into r leptons, we explored its indirect signals: one can reproduce the PAMELA e + 
excess and, at higher energies, a feature in the e + + e~ spectrum milder than the one hinted at 
by the ATIC data. 
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A Explicit realizations for Oqg and Oqf 
A.l Gluon operator 

To deconstruct Oqg one needs at least one heavy quark singlet F, of mass Mp and unspecified 
charge, and the mediation coupling XtpFF. Integrating it out at one loop gives rise to 

so that the threshold effect in Fig. [T]d could be obtained with M F ~ (A/47r) x 200 GeV. 
A. 2 Fermion operator 

As in the case of O5, there are two ways to generate Oqf, either by fermion or by heavy vector 
exchange. Concentrating on the vector-exchange case, the relevant gauge Lagrangian for a Z' of 
mass M could be 

L = g q Z'^ J2 nVi + 9^(<P + (A.2) 

i 

where i runs over all quarks. This can account for Oqf, with the needed effective scale to explain 
the tail of Fig.[T]D, provided 

M « (^/4vr) 1 / 2 v ^700 GeV. (A.3) 
At the same time, the Z'-width is dominated by decays into 00* : 



-,2 



T ~ —^-M. (A.4) 
48tt V ; 
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We believe that this leads to an acceptable situation, also in view of the dijet limits from CDF, 
provided — > Att. A model like this might be the easiest to incorporate in a complete extension 
of the SM. 
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